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ABSTRACT 


A  plethysmograph  was  developed  to  permit  non-invasive  measurement  of 
pulmonary  ventilation  in  small  animals.  In  this  system,  precise  measurements  of 
temperature,  humidity,  and  barometric  pressure  are  used  to  calculate  tidal 
volume,  while  a  BASIC  computer  program  directs  the  collection  and  display  of 
data  on  tidal  volume,  frequency,  body  temperature,  chamber  temperature,  and 
relative  humidity.  The  long  acoustical  time  constant  of  the  plethysmograph 
permits  accurate  measurement  of  tidal  volume  in  a  wide  variety  of  small 
unrestrained  animals  such  as  mice,  rats,  and  hamsters. 


INTRODUCTION 


The  two-chambered  plethysmograph  described  in  this  paper  permits  non- 
invasive  measurement  of  ventilatory  parameters  in  unrestrained  experimental 
animals  such  as  mice,  rats,  and  other  small  rodents.  Ventilatory  rate  and  tidal 
volume  are  derived  from  pressure  changes  within  the  animal  chamber.  Tidal 
volume  is  calculated  using  the  Drorbaugh  and  Fenn  equation  (5) .  This  equation 
corrects  for  temperature  and  vapor  pressure  changes  which  occur  as  inspired 
tidal  volume  is  warmed  and  humidified  and  as  expired  gas  is  cooled  and  loses 
water  vapor.  The  resulting  changes  in  pressure,  transmitted  throughout  the 
animal  chamber,  are  detected  by  sensitive  transducer  and  used  to  calculate  tidal 
volume . 

An  on-line  digital  computer  collects,  displays,  and  prints  out  all  data 
during  an  experiment.  The  software  program,  which  calculates  physiological 
parameters,  is  written  in  BASIC  and  was  developed  at  the  U.S.  Army  Research 
Institute  of  Environmental  Medicine.  The  program  which  has  automated  the  process 
of  collecting  data  from  large  numbers  of  animals,  collects  data  on  respiratory 
rate,  the  temperature  of  the  animal,  chamber  temperature,  relative  humidity,  and 
uses  this  information  to  calculate  tidal  volume  and  minute  ventilation. 

SYSTEM  DEVELOPMENT 
Description: 

The  plethysmograph  has  two  identical  chambers,  each  with  an  internal 
volume  of  11  liters.  Each  chamber  is  constructed  of  acrylic  plastic  and  has 
removable  0-ring  gasketed  lid  (Fig.  1).  Each  chamber’s  working  volume  is 
adjusted  to  10  liters  by  adding  1  liter  of  water.  The  addition  of  water  to  the 


chamber  insures  that  adequate  humidification  is  maintained.  Animals  are 
supported  above  the  water  by  a  plastic  grid.  Chamber  access  for  set-up  and 
cleaning  is  provided  by  a  removable  lid  which  is  secured  by  six-evenly  spaced 
90°  latch  clamps  which  provide  a  gas-tight  seal.  Temperature  and  humidity  are 
measured  by  a  dual-element  sensor  (#LIS-2052-PW,  Newport  Scientific  Co.).  A  set 
of  NBS  traceable  calibration  curves  are  supplied  with  each  sensor. 

A  balanced  ventilating  system  is  critical  for  the  operation  of  this 
system.  Ventilating  gases  are  introduced  into  the  chamber  via  a  pair  of  matched 
inlet  resistances  consisting  of  small  diameter  tubing.  The  resistances  are 
adjusted  to  provide  desired  flows  and  maintain  equal  pressures  in  both  chambers. 
The  inlet  resistances  eliminate  baseline  shifts  due  to  changes  in  temperature 
and  barometric  pressure.  The  desired  gas  composition  is  drawn  through  each 
chamber  by  a  regulated  vacuum  reservoir  (approx.  380  mmHg) .  The  outlet 
resistances  (Fig.  2,  section  #5)  are  adjusted  to  maintain  a  flow  of 
approximately  250  ml/min  and  equal  pressures  of  approximately  -1  cm  H2O  with 
respect  to  atmospheric  pressure  in  each  chamber  (Pbox) • 

For  calibration  purposes,  1  ml  of  ambient  air,  rapidly  injected  into  the 
box,  should  result  in  a  rise  in  Pbox  0.1  cm  H2O.  A  3.4  sec  acoustical  time 
constant  was  chosen  (Cgas  equals  lOml/cmf^O  and  Rinlet  0.34  cml^O/ml/sec) , 
since  it  is  satisfactory  for  the  measurement  of  respiration  in  rodents  and  other 
small  animals  with  high  ventilatory  rates  and  low  CO2  production. (1,2,4,5,13) . 

Layout : 

The  plethysmograph  is  composed  of  two  identical  acrylic  chambers:  a 
reference  and  an  animal  chamber  (Fig.  2).  The  2.54  cm  acrylic  wall  thickness 
prevents  distortion  of  the  chamber  (Fig.  2,  sections  #1  and  #2). 


The  pressure  differential  is  measured  with  a  transducer  (Validyne  MP45-1, 
range  *2  cm  H2O)  (Fig.  2,  section  #8)  and  amplified  with  an  amplifier  (Validyne 
CD-12) .  The  gain  is  adjusted  so  that  a  full-scale  pen  deflection  (50  mm  on  the 
strip  chart)  represents  a  pressure  change  of  0.125  cm  H2O.  Temperature  and 
humidity  measurements  are  displayed  by  a  digital  thermohygrometer.  Analog  D.C. 
signals  for  oxygen,  carbon  dioxide,  and  volume  are  displayed  on  a  recorder 
(Gould/Brush  Instruments,  Model  200).  A  mass  spectrometer  (Perkin  Elmer  Model 
1100A)  samples  gas  at  a  rate  of  60  ml/min  for  measurement  of  oxygen  and  carbon 
dioxide  fractions.  The  slow  leak  allows  a  gas  sample  to  be  taken  downstream  of 
the  system’s  last  resistor  (Fig.  2,  section  #10) .  The  coarse  pressure  balance 
in  the  box  is  maintained  by  an  adjustable  flow  meter , (Matheson  Model  7640  with 
Tube  #602)  (Fig.  2,  section  #9)  and  the  fine  pressure  balance  is  controlled  by  a 
micro-adjustable  valve  (Nupro  Model  B-4LA) (Fig.  2,  section  #3).  The  slow  leak 
of  gas  into  the  box  is  used  to  stabilize  box  pressure  (Fig. 2,  section  #10) . 

Solenoids  used  for  flushing  and  venting  can  be  activated  to  rapidly  change 
gas  composition  in  the  boxes.  Solenoids,  (Asco  Red-Top  TM?  1-1/4"  NPT) ,  are 
used  as  flush  valves  to  change  environmental  gas  mixtures  (Fig. 2,  section  #6) 
and  as  gas  input  valves  (Asco  Red-Top^  solenoid  valves,  3/8"  NPT) (Fig. 2, 
section  #7) .  Solenoid  and  gas  mixture  control  boxes  are  used  to  control  the 
sequence  of  events  and  composition  of  gases  in  the  chambers. 

Review  of  Barometric  Method: 

Plethysmography  requires  two  chambers,  one  for  the  animal  and  one  as  a 
reference  connected  by  a  differential  pressure  transducer  (1-5,7,8,11,12).  A 
slow  leak  mechanism  minimize  ambient  pressure  disturbances.  If  there  is  a 


change  in  ambient  pressure  the  slow  leak  mechanism  transmits  these  changes  with 
equal  magnitude  and  phase  shifts  to  each  chamber  without  effecting  the 
differential  pressure  measurements. 

Early  systems  were  unable  to  maintain  stable  baselines  due  the  inadequate 
chamber  sealing  (4,5).  Bartlett  and  Tenney  (2)  attempted  to  alleviate  unsteady 
baseline  problems  by  connecting  the  experimental  chamber  to  the  reference 
chamber  via  a  high  resistance  leak.  This  design  change  eliminated  some  baseline 
drift  found  in  closed  systems.  However,  earlier  continuously  purged  systems  had 
short  time  constants  which  resulted  in  a  significant  underestimation  of  the 
tidal  volume.  In  contrast,  Fappenheimer ’ s  design  (12)  has  high  resistance  inlet 
and  outlet  ports  and  requires  a  pressurized  gas  to  the  inlet  as  well  as  a  vacuum 
to  the  outlet  to  produce  a  flow  through  the  chamber.  This  and  other  similar 
designs  were  plagued  by  poor  frequency  responses,  excessive  drift  and  distorted 
respiratory  patterns  due  to  excessive  or  rapid  leakage  from  the  chamber. 

The  system  described  here  employs  a  design  which  is  very  similar  to  that 
of  Jacky  (8),  except  that  the  plethysmograph ’ s  flow  rate  is  greater  than  that  of 
Jacky’s  system  (8),  and  its  resistance  to  gas  flow  is  less  than  4  cm  H2O  with 
respect  to  atmosphere.  The  present  system  also  has  individual  resistance- 
inductance  networks  and  provides  smaller  inlet  resistances  and  a  greater  flow 
through  the  chamber.  A  regulated  vacuum  reservoir  allows  a  ventilating  flow 
through  the  chamber  which  is  adequate  for  use  with  300-600  gram  rats  who  produce 
a  moderate  amount  of  CO2.  This  level  of  flow  (250-270ml/min)  ensures  that  the 
CO2  concentration  does  not  rise  in  the  box.  This  system  differs  from  airtight 
chambers  because  it  behaves  as  a  first-order  high-pass  filter  with  a  corner 
frequency  of  5.3  cycles/min  and  has  a  flat  frequency  response  above  50 


cycles/minute.  The  3.4  sec  time  constant  is  large  enough  to  minimize  errors  and 
provide  an  accurate  representation  of  tidal  volume  through  the  range  of 
breathing  frequencies  exhibited  by  rodents.  The  reference  and  animal  chambers 
are  provided  with  identical  inlet  and  outlet  pathways  so  that  ambient  pressure 
disturbances,  some  which  are  similar  in  duration  to  those  of  the  respiratory 
cycle,  are  applied  with  equal  magnitude  and  phase  shifts  to  both  chambers. 
Therefore,  only  the  peak-to-peak  pressure  changes  which  occur  when  the  animal  is 
ventilating  are  observed. 

Calculations : 

Drorbaugh  and  Fenn  (5)  derived  the  following  formula  to  calculate  tidal 
volume : 

VT  =  Pt  .  Vk  .  Tr(Pb-Pc) 

Pk  Tr(Pb-Pc)  -  Tc(Pb-Pr)  where: 

Vf  =  tidal  volume  (ml) ; 

Vk  =  1  cc  volume  injected  into  chamber; 

Pt  =  pressure  deflection  associated  with  each  tidal  volume  (cm  H2O) ; 

Pk  =  pressure  deflection  associated  with  the  injection  of  calibration  volume  (cm 

H20); 

Tr  =  body  temperature  (rectal  temperature,  °C) ; 

Pc  =  vapor  pressure  of  water  in  the  chamber  at  room  temperature  (mmHg) ; 

Pb  =  barometric  pressure  (mmHg) ; 

Tc  =  air  temperature  in  animal  box  (°C) ;  and 

Pr  =  vapor  pressure  of  water  at  body  temperature  (mmHg) . 


There  are  several  significant  advantages  to  the  type  of  system  (12) 
described  in  this  paper:  1)  a  constant  flow  is  established  with  a  minimal 
pressure  drop  across  the  inlet  pathway;  2)  no  high  pressure  gas  source  or  fine- 
tuning  is  needed  to  produce  chamber  pressure  close  to  barometric  pressure;  3) 
there  is  no  noise  due  to  turbulent  gas  flow  caused  by  forcing  a  moderate  flow  of 
gas  through  inlet  orifices  under  pressure,  4)  the  acoustical  time  constant  of 
3.4  sec  is  long  enough  to  avoid  effecting  the  respiratory  pressure  signal  in  the 
range  of  frequencies  encountered  in  rats;  and  5)  the  modified  barometric  method 
(8)  allows  continuous  purging  of  the  chambers.  This  purging  reduces  CO2  levels 
to  0.02%  (PlC02)  under  steady  state  conditions  and  eliminates  errors  in 
ventilatory  measurements  caused  by  carbon  dioxide  build-up  in  the 
plethysmograph.  This  approach  avoids  the  problems  with  systems  which  required 
frequent  interruption  of  experiments  to  eliminate  carbon  dioxide  build-up  in  the 
chamber . 

Frequency  Response  Testing: 

The  barometric  method  has  been  compared  by  previous  investigators  (2,5,11) 
with  conventional  pneumotachographic  techniques  and  was  found  to  be  useful  as 
long  as  the  researcher  is  aware  of  known  artifacts  in  this  method.  Although 
most  investigators  did  not  purge  their  chambers  continuously,  this 
plethysmograph’ s  modifications,  which  include  the  gas-flush  component,  does  not 
affect  the  pressure-volume  relationship  of  the  plethysmograph.  Jacky  (8) 
determined  the  dynamic  response  of  his  plethysmograph  by  injecting  a  bolus 
volume  of  gas,  measuring  the  immediate  pressure  increase,  and  calculating  a  time 
constant  of  1.8  s  from  the  resulting  decay.  This  method  was  also  used  to 


measure  the  time  constant  (frequency  response  testing) (FRT)  in  this  system.  The 
excessively  short  acoustical  time  constant  in  Jacky’s  system  may  result  in  an 
underestimation  of  respiratory  volume  at  elevated  breathing  frequencies.  Thus, 
the  time  constant  in  the  our  system  was  lengthened  to  accommodate  awake 
unrestrained  animals  whose  breathing  frequencies  are  significantly  greater  than 
that  of  a  sleeping  mammal.  Our  calculation  of  the  tidal  volume  using  Jacky’s 
method  (FRT)  assured  that  our  system  also  had  a  flat  frequency  response  over  the 
entire  ranges  of  frequencies  encountered  in  rodents  under  a  variety  of 
conditions . 

Data  Acquisition  and  Analysis: 

Semi-automated  data  collection  and  analysis  is  performed  by  a  MINC  11/23 
computer  (Digital  Equipment  Corporation) .  Software  consists  of  three  chained 
programs  which  provide  calibration,  data  acquisition,  and  data  analysis 
functions,  respectively  (Appendix). 

The  first  program  allows  calibration  of  the  three  A/D  channels  for  volume 
change,  carbon  dioxide  and  oxygen  tension.  All  three  channels  can  be  calibrated 
by  providing  a  baseline  voltage  (volume  reference,  zero  for  PCO2,  and  PO2)  and 
then  injecting  a  known  calibration  volume  of  1  ml  to  measure  the  voltage  for  the 
volume  channel  as  well  as  known  gas  concentrations  to  be  used  as  gains  for  PCO2 
and  PO2.  Respective  calibration  factors  are  then  determined  and  stored. 

The  second  program  performs  the  data  acquisition  and  consists  of  three 
subroutines.  The  first  prompts  the  operator  for  a  data  file  name,  the  animal’s 
rectal  temperature,  barometric  pressure,  box  temperature,  and  relative  humidity. 
The  second  subroutine  provides  a  CRT  strip  chart  display  of  the  animal’s 


respiratory  pattern  (ie.  box  pressure)  and,  when  signaled  by  the  operator, 
collects  10  sec  of  data  at  a  digitizing  rate  of  50  Hz.  The  third  subroutine 
stores  the  data  collected  through  the  A/D  channels  and  keyboard  entries.  The 
program  allows  three  opportunities  for  data  collection  per  experimental  run. 

The  third  program,  which  analyzes  collected  data,  has  three  subroutines. 
After  calibration  file  and  data  file  identification  are  made,  volume  data  are 
graphically  displayed.  Ten  consecutive  breaths  can  be  selected  for  analysis. 
Tidal  volume  (Vf  (BTPS))  of  each  of  the  breaths  is  determined  using  the 
Drorbaugh  and  Fenn  equation  (5).  Mean  tidal  volume (V^),  respiratory  frequency 
(f ) ,  and  minute  ventilation  (Tf;)  are  calculated.  Mean  PCO2  and  PO2  are  also 
obtained.  A  hardcopy  printout  is  provided  of  all  identification  and  respiratory 
variables  measured  and  calculated. 

DISCUSSION 

Several  investigators  have  noted  limitations  in  previous  plethysmographic 
systems  (6,8).  These  have  included  an  unsteady  baseline,  fluctuations  in 
temperature  and  vapor  pressure  equilibrium,  excessive  ventilation  due  to  CO2 
accumulation,  and  inaccurate  measurements  of  differences  between  body  and 
ambient  temperatures. 

The  present  system  overcomes  most  of  these  limitations  by  providing  a  flow¬ 
through  design  with  individual  resistance-inductance  networks  compensating  for 
the  baseline  drift.  Because  the  respiratory  signal  is  so  small  it  is 
susceptible  to  distortion  caused  by  decreases  in  pressure  associated  with 
compression  of  intra-abdominal  gas.  This  problem  fortunately  produces  an 
easily  recognizable  elongation  of  the  expiratory  pressure  signal. 


Continuously  purged  plethysmographs  also  avoid  problems  with  temperature 
and  humidity  gradients  in  the  chamber  which  can  occur  in  studies  lasting  more 
than  10  minutes.  Epstein  and  Epstein  (6)  have  recently  shown  that  errors  in  the 
measurement  of  temperature  in  the  plethysmograph  are  responsible  for  even 
greater  underestimation  of  tidal  volume  and  developed  equations  for  correcting 
the  tidal  volume  calculation.  Tidal  volume  calculated  by  the  method  of 
Drorbaugh  and  Fenn  (5)  assumes  the  pressure  in  the  animal  chamber  increases  in 
proportion  to  V*p  as  inspired  air  is  warmed  and  humidified  during  inspiration. 
During  expiration,  pressure  decreases  with  cooling  and  water  condensation  from 
expired  gas.  Epstein  and  Epstein  (6)  suggest  that  the  pressure  changes 
associated  with  warming  and  humidification  of  inspired  air  and  the  subsequent 
cooling  and  water  condensation  are  not  equal.  Errors  of  15  to  30%  during  normal 
ventilation  and  10%  during  hypercapnic  exposures  are  reported  (6) . 

Our  calculations  are  in  agreement  with  Jacky  (9)  in  that  the  magnitudes  of 
the  errors  in  tidal  volume  are  less  than  those  of  Epstein  and  Epstein  (6)  who 
have  used  larger  animals  and  infants.  The  difference  in  measuring  tidal  volume 
by  conventional  pneumotachographic  and  plethysmographic  methods  have  been  the 
subject  of  much  research  (6,7,9,10)  with  each  group  presenting  different 
solutions  for  correcting  the  error  in  the  calculation  of  tidal  volume. 
Plethysmographic  measurements  of  tidal  volume  have  been  criticized  by  Epstein 
and  Epstein  as  having  a  greater  underestimation  of  true  tidal  volume  than 
conventional  methods  such  as  pneumotachographic  or  mechanical  methods 
(2,4,5,11).  They  believe  that  the  asymmetry  of  the  pressure  events  (of 
inspiration  and  expiration  and  associated  time  constants)  leads  to  an 
overestimation  in  inspiratory  time  (Ti)  as  the  total  breath  duration  (Ttot) 


increases.  Epstein  and  Epstein  contend  that  the  assumptions  of  other 
investigators  that  expiration  mirrors  inspiration  is  incorrect.  However, 
Epstein’s  finding  are  not  supported  by  the  literature  (2,4,5,11)  and  we  believe 
that  the  errors  which  are  inherent  in  other  plethysmographs  are  eliminated  by  a 
"leaky"  system.  When  comparing  resting  data  from  a  recent  study  with  rats  with 
that  of  stimulated  breathing  (13)  we  found  no  large  errors  in  ventilation. 
Jack y  (9) ,  has  reported  that  the  approximate  magnitude  of  potential  error  in 
tidal  volume  measurement  is  in  the  order  of  15%  underestimation  during  air 
breathing  and  decreases  to  4%  underestimation  of  the  fractional  increase  in 
tidal  volume  during  hyperpnea  due  to  hypercarbia,  associated  with  changes  in 
timing  of  the  respiratory  cycle  and  is  varied  with  different  breathing  patterns 
(10,13).  Our  measurements  were  not  corrected  for  the  underestimation  in  tidal 
volume  because  these  errors  are  only  significant  when  absolute  measurements  of 
tidal  volume  are  made,  but  not  when  relative  changes  were  considered  as  in  our 
study  (13) .  Typically  our  acute  studies  did  not  require  using  any  one  gas  for 
longer  than  5-7  minutes  and  this  duration  is  short  enough  to  minimize 
temperature  distortion  of  the  respiratory  signal. 

Practical  and  theoretical  considerations  must  be  reviewed  when  the 
plethysmograph  is  used  in  research  studies.  The  plethysmograph  must  be  well 
constructed  and  undergo  rigorous  testing.  The  investigator  must  be  familiar  with 
the  system  because  the  small  respiratory  signal  is  easily  distorted  by  internal 
and  external  factors.  Once  these  considerations  are  addressed,  the  operation  of 
the  system  is  suitable  for  measuring  ventilation  in  awake , unrestrained  animals. 
The  following  are  critical  characteristics  and  components  of  the 


plethysmographic  system:  1)  the  external  reference  chamber,  2)  a  stable  and 


vibration-free  base,  3)  a  relatively  long,  time  constant  to  damp  out  slow  changes 
in  pressure  due  to  temperature  and  humidity  effects,  4)  CO2  scrubbers  or  a 
continuously  flushed  "leaky"  system  to  avoid  CO2  build  up,  and  5)  an  effective 
system  for  measuring  the  partial  pressures  of  carbon  dioxide  and  oxygen. 
Although  more  direct  methods,  such  as  pneumotachography ,  have  been  shown  to  have 
smaller  errors  in  measuring  "true"  tidal  volume,  the  plethysmograph  is  clearly 
suitable  for  monitoring  changes  in  ventilation  in  small  unrestrained  animals. 


SUMMARY 


A  plethysmograph  system,  with  an  automated  gas-control  system  has  been 
developed  for  the  non-invasive  measurement  of  pulmonary  ventilation  in  small 
animals.  A  relatively  long  time  constant  preserves  signal  fidelity  and  removes 
drift  while  a  slow-leak  system  limits  the  effects  of  ambient  pressure 
disturbances.  The  proper  selection  of  the  time  constant  and  accurate  measurement 
of  temperature  and  humidity  reduces  the  underestimation  in  tidal  volume  in  the 
plethysmograph.  The  ease  with  which  the  gas  mixtures  can  be  administered  is  also 
conducive  to  longer  and  more  variable  study  designs.  Early  problems  with  carbon 
dioxide  buildup  are  easily  eliminated  by  using  a  continuously  flushed  system 
whose  ventilation  flow  rate  is  set  at  250-270  ml/min.  Precise  measurement  of 
temperature  and  humidity  in  the  animal  chamber  permit  accurate  tidal  volume 
corrections.  In  conclusion,  the  plethysmograph  design  presented  in  this  paper 
reflects  several  years  of  design  and  research  to  develop  a  functional  tool  to 
measure  pulmonary  ventilation  in  unrestrained  small  animals. 
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APPENDIX 


Computer  Program: 


£8  Bffi  x ( 2  > 

30  DI S PL A V _CLEAR 
<0  PRINT  'OFT:  A 
50  PRINT  'tree! 

60  PRINT  ' 

70  PRINT  ' 

60  PRINT  ' 

90  IF  QX>  =  1  THEN 
100  C.0  TO  30 
110  ON  Q7.  GO  TO  1 
120  R£t1*******«»* 


PROGRAM  I  VFT************************* 

PROGRAM  FOR  DATA  ACQUIST10N  &  ANALYSIS' 

1)  calibrate!' 

2)  collect  data!' 

3  )  ana  lyre  data ! ' 

4)  *  top  !  '  #  \  I  NPUT  ij'/. 

IF  U*t=4  GO  TO  !!0 

30. ! 40 . 1 50 • 160 

********l i st  of  chained  programs******* 


\  PRINT 


\  PRINT 


130  GO  TO  170 
140  CRAIN  ' OF  1 T  ' 

150  CHAIN  ' OF27 ' 

160  GO  TO  1220 

170  REM*******«****«***CAL1BRATI0N  Or  VOLUME*************®* ******** 

180  DISPLAY-CLEAR 

190  PRINT  'A/D  0  VOLUME' 

200  PRINT  A/D  1  02' 

210  PRINT  'A/D  2  C02 '  \  PRINT  \  PRINT 

220  PRINT  'CALIBRATION  OF  VOLUME'  \  PRINT 

230  PRINT  'type  1)  to  enter  calibration  factor;' 

240  PRINT  '  2)  to  calculate  calibration  Factor:':  \  INPUT  07. 

250  IF  0X>*1  THEN  IF  QX<=2  GO  TO  270 

260  C.0  TO  170 

370  IF  dX*2  GO  TO  300 

290  PRINT  \  PRINT  '  enter  calibration  factor!';  \  INPUT  F3 
290  GO  TO  650 
300  DISPLAY-CLEAR 
310  DIM  V< 100) 

320  PRINT  'DETERMINATION  OF  CALIBRATION  FACTOR  BY  INJECTION  OF  1  ML' 
330  FOR  1*1  TO  5 

340  PRINT  \  PRINT  'CiiECK  BASELINE'  \  PRINT 
350  PRINT  'when  reaoy  press  RETURN'! 

360  INPUT  X* 

370  V*0 

380  AIN< .VI ) .40. .05) 

330  FOR  11=0  TO  33 


400  V*V*V( I  1  ) 

410  NEXT  II 
420  BCD  =V/40 

430  PRINT  'INJECT  1  ML  NOW ' ' ! CHRS C 7 ) 

440  A I N ( , V( ) . 100. .05) 

450  P*B ( I )*1 .02 
460  FOR  12=0  TO  99 
470  PRINT  12. V< 12) 

4eO  IF  ABS  <  V ( I  2 )  ) >  =  P  THEN  P  =  A9S < V < I  2 ) > 

490  NEXT  12 
500  F  < I ) *P-e< I  ) 

510  PRINT  'BASELINE  ( uo 1 1 s  )  I  '  !  B  (  I  )  .  '  CAL  #  '  >  I  !  '  ( uo 1 1 s /m 1 ) I  ' I F (  I  ) 
520  NEXT  I 

530  REM* »**••*•»» AVER AGE  OF  CALIBRATION  INJECTIONS********** 

->40  F  3  =  0 

550  PRINT  'do  you  wish  to  delete  a  calibration  (Y/N)'»  \  INPUT  G* 
560  IF  US  =  ' N  GO  TO  530 

570  PRINT  'hit  number  of  run  to  delete'!  V  INPUT  RS 
590  13=0 

»30  FOR  1*1  TO  5 
600  IF  I =  R9  THEN  630 
S’.O  F3-F3+-  (  I  ) 

520  13=13+1 
630  NEXT  I 
640  F3  =  F  3/ I  3 

650  PRINT  \  PRIM  'CALIBRATION  FACTOw  <VOlTS/ML>  *  ' ! r 3 
360  PRINT  \  PRINT  'do  you  with  to  recalibrate  (Y/N)';  \  1NDUT  0* 
S  70  IF  0*=  Y'  GO  TO  170 
580  X ( 0 ) *F 3 

680  DI6PCRY  CLEAR 
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710  IF  H=  1 0  GO  TO  7S0 

720  FOR  K*J  TO  S2X 

720  IF  BIKKP  THEN  P*B(K) 

740  IF  B ( K ) >P  GO  TO  760 
730  NEXT  K 
760  P0  =  P  \  I  1 =  K 
770  P-B<  1 1  > 

7B0  NEXT  II 

790  FOR  12=1  TO  10 

BOO  0=0 +  V( 12) 

810  NEXT  12 
820  0=0/10 

830  F  =600/  (  (  S27.-S  1  '/.)*.  02  ) 

840  FOR  11=511  TO  609  STEP  2 
850  V6=V6+A3S(B<  ID) 

860  V7=V7+A9S(B(  I  1  +  1  )  ) 

870  NEXT  II 

880  06  =  ABS( (06/50) /X<  1)  1 
890  07=ABS< (07/50)/X<2) ) 

900  DISPLAY.CLEAR 

910  PRINT  'RESPIRATORY  RATE  =  'IF 

920  PRINT  'TIDAL  VOLUME  =  ' ! V 

930  PRINT  'MINUTE  OENT I L AT  1  ON  =  DF*V 

940  PRINT  ' P02  =  '  106 

950  PRINT  ' PC02  =  '107 

960  GOSUB  1010 

970  PRINT  \  PRINT 

980  PRINT  'to  continue  data  analysis  hit  RETURN'! 
990  NEXT  I 
1000  GO  TO  140 

1010  REM»***«**«*PR1NTF.R  SUBROUTINE*************** 
1020  OPEN  'IP '  FOR  OUTPUT  AS  FILE  #1 
1020  PRINT  #1  \  PRINT  41 

1040  PRINT  41,DAT*,CLK* 

,1050  PRINT  41, 'FILE  ID  =  SYlI'FtiN* 

1060  PRINT  41. 'RAT  NO.  =  DA* 

1070  IF  A 1  *  =  '1  '  THEN  PRINT  41,  'NO  DRUG  -  N0R10X1A' 


\  INPUT  D* 


THEN  PRINT  41 
THEN  PRINT  #1 
THEN  PRINT  41 
THEN  PRINT  #1 
RUN  NUMBER  «’  ' 


THEN  PRINT  #1 
THEN  FRINT  41 
THEN  PRINT  41  , 


1080  IF  A!*='2 
1020  IF  A 1  *  =  1  3 
1100  IF  A 1 *=  '4 
1110  IF  Al*  = '5 
1120  IF  A 1  *  * '6 
1130  IF  A 1 ♦ = '7 
1140  IF  A 1 1  =  '8 
1150  PRINT  41  , 

1160  IF  A3*  =  ' 0 
1170  IF  A3*=  '  1 
1180  IF  A3*=  '2 
1190  IF  A3*  =  '2 
1200  IF  A3*  =  '4 
1210  PRINT  #1, 

1220  PRINT  41 , 

1230  PRINT  41 , 

1240  PRINT  41 , 

1250  PRINT  41 , 

1260  PRINT  41 , 

1270  PRIM  41  , 

1260  PRINT  41, 

1290  PRINT  41 , 

1200  PRINT  41 , 

1310  PRINT  41 . 

1320  PRINT  41 
1230  FOR  : 2- 1  IQ  lo 
1240  PRINT  41,0(13) 

1250  NEXT  13 

1260  PRINT  41  \  PRINT  41 

1370  PRINT  4) ,X(o> ,X( 1 ) ,X(2> 

1360  CLOSE  41 

1390  RETURN 

1400  DISPLAY. CLEAR 

14  10  '.''.OSE 

\\Ub  U40 


THEN  PRINT  41 » 'NO  DRUG  -  HYPOXIA' 

THEN  PRINT  41,  'NO  DRUG  -  HYPOXIC,  HYPERCAPNIA  ' 
THEN  PRINT  41, 'NO  DRuG  -  HYPERCAPNIA' 

'DRUG  -  HYPEROX  I C  HYPERCAPNIA ' 

-  HYPOXIA' 

-  NORi+OX  I A  ' 

-  HYPERCAPNIA' 


.  'DRUG 
,  'DRUG 
,  'DRUG 
:  A2* 


'NO  DRUG' 

'DOSE  =  0.005  ML/KG ' 

'DOSE  =  0.050  ML/KG' 

'  THEN  PRINT  41,'D0S€  =  0.500  ML/KG' 

'  THEN  PRINT  41, 'DOSE  =  5.000  ML/KG' 

'RECTAL  TEMP  ( C )  =  '  ! B ( 1 > 

'eOX  TEMP ( C )  =  ' B ( 2 ) 

'BAROMETRIC  PRESSURE  <  TORR )  =  'BO) 
'RELATIVE  HUMIDITY  *  ’B(4) 

'PH20  (BOX)  =  ' E ( 5 ) 

'  PH20  (RAT)  =  ' ! B ( 6 ) 

'TIDAL  VOLUME  =  '!V 

'RESPIRATORY  RATE  =  '  F 

'MINUTE  VENTILAI10N  = 

'  rU2  = 


' PCQ2  = 
\  PRINT 


I  F  »V 
• :  vs 
■ :  v7 


41 


j 

J 

‘4 

* 


| 

i 
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N 

? 

,\ 
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10  REM* ********** PROGRAM  NAME  -VF1T-********* 
20  REM***CHANNEL  0  PRESSURE  ,  CHI  02.  CH2  C02 
30  DIM  M(SOO) .G< 100) .8(610) 

*0  DISPLAY-CLEAR 


30  PRINT  'RAT  NO.  »'! 
60  PRiNT  'CONDITION 
70  PRINT  ' 

80  PRINT  ' 

90  PRINT  ' 

100  PRINT  ' 

110  PRINT  ' 

120  PRINT  ' 

130  PRINT  ' 

140  PRINT  ' 

150  INPUT  Al* 

160  PRINT  'RUN  #  (0-S) 
170  PRINT  'DOSE 


\  INPUT  A* 

1-NQ  DRUG  -  NURHOXIA' 
2* NO  DRUG  -  HYPOXIA' 


2* NO  DRUG  -  HYPOXIA' 

3- NO  DRUG  -  HYPOXIC.  HYPERCAPNIA' 

4 - NO  DRUG  -  HYPERCAPNIA' 

S-NO  DRUG  -  HYPER0X1C  HYPERCAPNIA' 


6*DRUG 

7- DRUG 

8- DRUG 

9- QTHER ‘ 


-  HYPOXIA' 

-  NORMOXIA' 

-  HYPERCAPNIA' 


■I  \  INPUT  A2* 

170  PRINT  'DOSE  0*  NO  DRUG' 

180  PRINT  '  1  * . 005  ML/KG' 

190  PRINT  '  2- . 050  ML/KG' 

200  PRINT  '  3*. 500  MI./KG' 

210  PRINT  '  4-5.00  ML/KG' 

220  INPUT  A3* 

230  Al*-  'SY1 I  '&AS&A1 S&A2S&A3S 
240  PRINT  'DATA  SAVED  ON  FILEI'Al* 

250  PRINT  'RECTAL  TEMP  (C)  «  ' I  \  INPUT  B( 

260  PRINT  'BOX  TEMP  (C)  -  '!  \  INPUT  BC 

270  PRINT  'BAROMETRIC  PRESSURE  < TORR )  -  ' J  \  INPUT  BC 

280  PRINT  'RELATIVE  HUMIDITY  »  'I  \  INPUT  BC 

290  REM  CALCULATION  OF  PH20 

300  C* 1 . 23234E+07  \  Cl-7. 52247  \  C2— 1223.31  S  C3--22 
310  C6*C1  +  (C2/C4)  +  (C3/(C4'>2)  ) ♦ ( C/ (C4A3 )  )  \  BtSl-lO'CS 
320  C8  =  C1*(C2/C5)-MC3/(C5A2))  +  (C/(C5'>3>>  \  B(6)-10ACB 
330  B(7) *0  \  B ( 8 ) -0  \  B(9)=0 


INPUT  B ( 1 ) 
INPUT  B ( 2 ) 
INPUT  B ( 3 > 
INPUT  B ( 4 ) 


C2— 1223.31  S  C3--222614  \  C4-B(2>+273  \  C5-B(l)+273 


340  REM 

350  DISPLAY- CLEAR 
360  PRINT  'typ#: 
370  PRINT  ' 

380  PRINT  ' 


1)  new  File  information' 

2)  collect  data' 

31  return  to  MAIN  MENU ' 1 


\  INPUT  G7. 


390  IF  GX>» 1  THEN  IF  QX<*3  GO  TO  410 
400  GO  TO  350 


410  ON  OX  GO  TO  40,420.760 

420  REM*******»»**DATA  COLLECTION  ROUTINE******************** 

430  FOR  1-1  TO  3 

440  IF  1*1  THEN  N*« '1 ' 

450  IF  1-2  THEN  N«='2' 

460  IF  1-3  THEN  N*-'3' 

470  OPEN  Al*»N*i'.DVF'  FOR  OUTPUT  AS  FILE  #2 
480  AIN<  'DISPLAY' .M< ), 1 ,  ,0) 

490  PRINT  'DATA  FILE  '  Al *&N* 

500  PRINT  'press  RETURN  when  readr  to  sample.' 

510  AIN(  'DISPLAY  '  ,M( ). 1 , ,0) 

520  GET-CHAR (C*l 
530  IF  LEN ( C* ) -0  GO  TO  510 
540  IF  ASC <  CS  )  =  1 0  GO  TO  560 
550  GO  TO  510 

560  AIN( ,M( ) .500, .02)  \  A IN ( . G ( ) , 1 00 , . 02 . 1 . 2 )  \  REM  SAMPLE  CHANNELS  0-2 


570  PRINT  CHR*(7>! 

580  FOR  1 1-1  TO  3 
590  PRINT  #2 . B ( 1 1 ) 

600  NEXT  II 

610  FOR  11=10  TO  510 

620  B(I1)*M(I1-10) 

630  PRINT  #2 . B ( 1 1 ) 

640  NEXT  II 
650  FOR  11-511  TO  610 
660  B( II >*G< 1 1-51 1 ) 
670  PRINT  PZ.B(Il) 

680  NEXT  11 

690  CLOSE 

700  NEXT  I 

710  DISPLAY. Cl  EAR 

720  PRiNT  '  t  y p* l  1  ) 

730  PRINT  '  2) 

740  PRINT  '  J ) 


1)  new  file  information’ 

2)  return  to  MAIN  MENU' 

3)  stop-;  \  INPUT  87. 


750  IF  QX>*1  THEN  IF  Q'/.<-3  GO  TO  770 
760  GO  TO  710 


770  ON  OX  GO  TO  40.780.7S0 
780  CHAIN  ' VF  T ' 

790  END 


rag 


I  4'  M  |.<  <.|  1.1 


IFitTWit 


10  REM***********PRaGftAM  NAME  “UF2T -**«***«***# »*♦*»*•••*• 

20  REM***CHANN£L  0  PRESSURE .  CHI  02.  CH2  C02  •••******•***«* 

30  DIM  6(610  > ,X<2> . V( 1 1 ) 

40  D I  SPLAY_CL£AR 

30  PRINT  \  PRINT  \  PRINT 

60  REM*********0P£N  CALIBRATION  FILE**********»***«»* 

70  PRINT  'ENTER  CALIBRATION  FILE  NAME':  \  INPUT  C* 

80  OPEN  'SY1 :  '4CS4  '  .CAL  '  FOR  INPUT  AS  FILE  #2 
90  FOR  1=0  TO  2 
100  INPUT  #2 . X < I ) 

110  NEXT  I 

120  CLOSE  #2 

130  PRINT  \  PRINT 

140  PRINT  'CALIBRATION  FILE  ISI'IC* 

150  REM*********HAIN  MENU****** ******* ******* 

160  PRINT  \  PRINT  \  PRINT 

170  PRINT  'type!  1)  for  new  calibration  File' 

180  PRINT  '  2)  to  analyze  neu  data  File' 

190  PRINT.  '.  3)  tJB  .collect  neu  data' 

200  PRINT  '  4)  to  stop':  \  INPUT  Q'/. 

210  IF  GX>« 1  THEN  IF  0XY=4  GO  TO  230 
220  GO  TO  150 

230  ON  Q7.  GO  TO  50.240,250.1400 
240  GO  TO  260 
250  CHAIN  'OFT' 

260  REM** ****** **OPEN  DATA  FILES**********»**«* 

270  DISP^AY.CLEAR 
2B0  PRINT  \  PRINT 

2S0  PRINT  'ENTER  DATA  FILE  NAME':  \  INPUT  FS 

300  FOR  1=1  TO  3 

310  IF  1  =  1  THEN  N*= ' 1  ' 

320  IF  1=2  THEN  NS  =  '2' 

330  IF  1=3  THEN  NS  = '3' 

340  OPEN  'SY1 :  '4FS4NS4 ' .DUF  '  FOR  INPUT  AS  FILE  #3 
350  FOR  11=1  TO  610 
360  INPUT  #3.6(11) 

370  NEXT  II 

3eo  CLOSE  *3 

390  A*=SEG*(F*. 1,1) 

400  A1»=SEG*(F*.2,2) 

410  A2*=SEG*(FS.3,3> 

420  A3*=SEG*(F*.4,4) 

430  REM**********»**DATA  ANAL YS I S** »**»**•******«»»• ******* 

440  DISPLAY.CLEAR 
450  S 1 X  =  0  \  S2'/.  =  0 
460  GRAPH! ,500. ,e< 10) > 

470  LABEL ( »  'RAT  I  '4F*4iS) 

480  PRINT  'shaae  graph  (Y/N)','  \  INPUT  Q* 

490  IF  H*=  'N  '  C.0  TO  510 
500  SHADE 

510  PRINT  'select  beginning  cY  First  breath  (  -:  S 1  ’/.:')  1'  :  \  INJuT  Qi 
520  IF  l.£N(Q*K>0  THEN  31  X*UA._  I  3*  ) 

530  PRINT  'select  end  of  ten-.n  breath  (  '  I  SEX :  '  >  1  '  :  \  INPUT  OS 
540  IF  LENUISIOO  THEN  S2X  =  o«l  l  OS  ) 

550  PQ I  NT  (  '  ERAND  '  ,  SIX .  t  (  S IX. )  ) 

5S0  POINT (  'BRAND ' .  S  2  X  >  P  (  S  2  X  > ) 

570  PRINT  'regranh  points  (Y/N)':  \  IN^UT  OS 
Seo  IF  0 S  =  ' Y '  GO  TO  460 

590  PR  INI  'change  start  or  step  pcmt7  (Y/N)-:  \  INPUT  OS 
600  IF  3 S  =  '  Y  lib  TO  510 

610  RE -»«»***»»»***» -EAbUREMcNT  0-  TIDAL  OUbU'-E** »«*••*••»**« »»« 

620  PO=B (SIX)  \  P=e(S:X)  \  M  =  u  \  s>  =  0  \  F=0  \  06-0  \  07=0 

630  FOR  Il=s:x  TO  S2X  j 

640  IF  B< I  1 >  >  =  P  GO  TO  770 
650  M=M+1 
660  V(M)*P-PO 

670  0(M)=(0(ri)/X(0) )*( <B( 1 )+273)*(B(3)-B(5) ) )/< (B< 1 )*273)*(B(3)-B(5) )-(B(2)*27 
680  IF  O(M) >.5  GO  TO  700 
690  n=n-i 
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700  REM*e**e*#********#»CALIBRAT ION  OF  GASES**«*»«e»<************e 

710  PRINT  'CALIBRATION  OF  02  AND  C02 '  \  PRINT  \  PRINT 
720  PRINT  'type  I)  to  emir  calibration  factor!' 

730  PRINT  '  2)  to  calculate  calibration  facto  /  'I  \  INPUT  Q7. 


740  IF  (3%>*1  THEN  IF  G%<*2  GO  TO  760 

750  GO  TO  690 

760  IF  Q-/.»2  GO  TO  820 

770  PRINT  \  PRINT  'enter  calibration  factor  for  02  ( uo 1 ts/ tcrr )J  '  ! 
780  INPUT  X ( 1 ) 

790  PRINT  'enter  calibration  factor  for  C02  ( #c 1 ts/ torr >  I  ' ! 

800  INPUT  X <  2 ) 

810  GO  TO  1040 
820  DISPLAY.CLEAR 

830  PRINT  'DETERMINATION  OF  GAS  CALIBRATION  FACTORS' 

840  DIM  G<  20 )  \  N  =  0  \  G*0 

850  PRINT  \  PRINT  'sample  zero  02.  when  ready  press  RETURN'! 

860  INPUT  X* 

870  G0SU3  1150 

880  PRINT  'zero  02*  ' ! G 1 !  ' uo i t s  ' 

890  PRINT  \  PRINT  'sample  zero  C02,  when  ready  press  RETURN'! 

900  INPUT  X* 

910  GOSUE  1150 

920  PRINT  'zero  CU2*  '  ! G2 ! 'uo  1 1 s  ' 

930  PRINT  \  PRINT  'samrie  peaK  02.  when  ready  press  RETURN'! 

940  INPUT  X* 

950  G0SU3  1150 

960  PRINT  'peart  02*  '  !G3! 'uolts  ' 

870  PRINT  'enter  peaK  02  torr!'!  \  INPUT  D5 

980  PRINT  \  PRINT  'sample  peaK  C02.  when  ready  press  RETURN'! 

990  INPUT  X* 

1000  GOSUB  1150 

1010  PRINT  'peaK  C02*  '  !G4!  'uolts  ' 

.1020  PRINT  'enter  peaK  C02  torr!'!  \  INPUT  D6 
1030  X< 11=ABS< (G3-G1 ) /D5)  \  X ( 2 ! *A8S ( < G4-G2 ) /D6 ) 

1040  PRINT  \  PRINT  '02  CALIBRATION  FACTOR  ( UOLTS /TORR >*'! X I  1 ) 

1050  PRINT  \  PRINT  'C02  CALIBRATION  FACTOR  ( UOLTS/TORR ) *  ' ! X ( 2 ) 

1060  PRINT  \  PRINT  'do  you  wish  to  recalibrate  (Y/N)'!  \  INPUT  0* 
1070  IF  G*«  ' Y '  GO  TO  690 

,  1 080  PRINT  'ENTER  CAL  FILE  NAME  < da y . mon t h  ,  n umb er >  ' !  \  INPUT  N« 
■1090  G°EN  'SY1 I  '4N»4  '  .CAL '  FOR  OUTPUT  AS  FILE  #2 
1100  FOR  1*0  TO  2 
1110  PRINT  #2 • X (  I  ) 

1120  NEXT  I 
1130  CLOSE  »2 
J 1 40  GO  TO  30 
1150  N  =  N+1 

1160  IF  N* 1  THEN  D*1 

1170  IF  N*2  THEN  D=2 

1180  IF  N=3  THEN  D* 1 

1190  IF  N=4  THEN  D=2 

1200  AIN( .G< ) .20, . 1 ,D> 

1210  FOR  1*0  TO  19 
1220  G  =  G-G< I  ) 

! 1230  NEXT  I 

1240  IF  N= 1  THEN  G!=G/20 
1250  IF  N  =  2  THEN  G2*Ci/20 
1260  IF  N  =  3  THEN  G3  =  G/20 
1270  IF  N*4  THEN  lt4  =  G/20 
12H0  FUR  1*0  TO  20 
1290  It  (  I  I  * O 
1300  NEXT  I 
131"  RETURN 
1320  END 
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